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1 Introduction 

While inflation has become a very successful paradigm for early-universe physics, it still leaves 
us with many challenges and puzzles. Perhaps the most important of these is incorporating 
it into a broader framework of quantum gravity, which allows systematic treatment of initial 
conditions, quantum corrections, etc. Even lacking such a treatment, quite a bit of detailed 
study can be accomplished, treating general relativity as an effective field theory and analyzing 
it perturbatively. One of course faces the usual problem of nonrenormalizable UV behavior, 
but handles this by imagining that appropriate counterterms are provided by a more complete 
theory, and do not appreciably affect the long distance behavior of the theory relevant for 
inflationary predictions. 

However, another puzzle arises in the infrared behavior of the theory: various quantities 
in de Sitter space, or other inflationary cosmologies, exhibit IR divergences. IR divergences 
typically indicate that one has not understood an important aspect of the dynamics, or has 
asked a bad question. For this reason these have been extensively discussed, with arguments 
ranging from the statement that they can be innocuously eliminated to the statement that they 
indicate a fundamental instability of de Sitter space causing the cosmological constant to relax 
to zero [iHTJ. It is important to more clearly understand where such IR divergences occur, and 
in what way they modify calculation of physical quantities. 

While IR divergent effects of massless scalars have long been known (see e.g. [8|), one might 
question how general these are, given that the correct theory may not have massless scalars. 
However, in an effective description of gravity masslessness is inherent in a description in terms 
of gravitons. This strongly motivates understanding the IR behavior of processes involving 
gravitational loops. 

Loop calculations in gravity are notoriously complicated, due to the higher spin and its 
corresponding gauge invariance. For this reason, so far results have been very limited. One 
of the goals of this work is to improve this situation. In particular, we calculate gravitational 
loop corrections to propagators in de Sitter space. While we don't exhibit all such corrections, 
we note that some of the corrections are indeed IR divergent, and we extract this leading 
IR behavior. Moreover, we argue that there is a simple way to derive such behavior from 
semiclassical arguments, related to the consistency relations of pjjll] (and with antecedents 
in 



12 l3j). We explicitly check examples where these semiclassical relations match the detailed 
in-in loop calculations. 

We then propose extending the semiclassical methods to more general, slow-roll, cosmolo- 
gies, and show how they give leading IR corrections there. We also demonstrate another check 
of these methods based on an explicit calculation. Interestingly, we find that the corrections 
implied by these methods shift various cosmological observables, notably the tensor /scalar ra- 
tio and nongaussianity parameters in the "large box" (total inflated volume), and we exhibit 
leading contributions to such shifts. Since the corrections grow in the IR, for a sufficiently long 
period of inflation, for example in models with slow-roll after exit from a chaotic phase, the 
loop corrections become of order one, compared to the tree-level quantities. This suggests a 
breakdown of the perturbation series, associated with the strong IR effects. 



This apparently sharpens the question of the physical interpretation of IR growth in our 
descriptions of inflationary spacetimes. While there have been various suggestions that IR 
effects in de Sitter can be resummed, it seems possible that the breakdown of the perturbation 
series represents the need for new non-trivial non-perturbative dynamics. 

Indeed, another motivation for this work is an interesting parallel, in black hole physics. 
There are good indications that the extreme UV (superplanckian) regime of gravity is governed 
by physics classically approximated as black hole formation, providing a link between UV and 



IR physics. A critical question, in light of Hawking's perturbative arguments 14 for informa- 
tion loss, is how such physics can be unitary. Specifically, by a calculation of the state on a slice 
spanning the exterior and interior of a semiclassical black hole spacetime, one constructs an 



argument for information loss. However, arguments were given in 15 , 16 that the perturbative 
state on such a slice ceases to be a good description, by a sufficiently long time, t ~ RS, where 
R and S are the radius and the entropy of the black hole; one such argument was based on 
apparently large fluctuations in the perturbative state. Failure of such a sharp argument for 
information loss was proposed as a resolution of the information paradox, leaving an informa- 
tion problem of determining the non-perturbative gravitational mechanics that completes the 
description in this regime. Quantization on such slices in black hole spacetimes is very similar 
to quantization on spatial slices in inflationary spacetimes. The calculations of this paper ex- 
plicitly show perturbative corrections to the latter that become large on the time scale t ~ RS, 
where now R and S are the de Sitter radius and entropy, in close parallel to the black hole 
storyjj This deepens the analogy between the situations, and amplifies the arguments for a 
similar breakdown in the black hole context. 

In outline, the next section summarizes in the de Sitter (dS) case the basic semiclassical 
relations and resulting leading IR effects on propagators, due to tensor loops. Section three 
then summarizes checks of these results, through explicit loop calculations of in-in correlators, 
in which the leading IR behavior is extracted. Section four extends the semiclassical relations to 
slow-roll spacetimes, accounting for both tensor and scalar loops, and performs another check 
on the formalism based on results in the literature. There is also a brief comment on the relation 
to the 6N approach. Section five calculates the leading corrections to cosmological observables, 
such as the tensor to scalar ratio, and nongaussianity parameters. It also contains a worked 
example of scenarios of slow-roll following chaotic infiation, where such corrections become large 
due to long-time evolution. Section six closes with further discussion of interpretation of these 
IR effects, and the parallels with black hole physics. The appendices contain review of previously 
known consistency relations, more complete derivation of the semiclassical consistency relations 
that we use, and more details of the in-in loop calculations. 



^This analogy was in particular pursued in ref. jlTJIlS], where it was pointed out than an attempt to regulate 
inflation by having it terminate through slow-roll fails if the termination time exceeds t ^ RS. The mechanism 
for this failure is the transition to self-reproduction. 



2 Semiclassical relations and IR effects in de Sitter 

We first consider correlators and their loop corrections in the simple case of de Sitter space, with 
a free scalar field a of vanishing expectation value coupled to the metric. Here the background 
is fully determined by the cosmological constant A. The action is 

S=^JV^[R- d.ad'-a - 2A] , (2.1) 

where R is the spacetime Ricci scalar, and we choose units 8ttG = l/iW? = 1. A perturbative 
description of the coupling between matter and metric can be derived using the Arnowitt- 



Deser-Misner (ADM) formalism 19 , which gives the action for the coupled perturbations to 



any given order by an iterative procedure. The ADM form of the line element is 

ds"^ = -N'^dt^ + hij{dx^ + N'dt){dx^ + N^dt) , (2.2) 

where iV, N'^ are the lapse and the shift functions. The three-dimensional metric hij encodes 



scalar and tensor fluctuations in the spatial geometry. When the line element (2.2) is inserted 



into Eq. (2.1), the lapse and the shift functions act as Lagrange multipliers: the field equa- 
tions obtained by extremizing the action in N and A^' give the constraint part of the Einstein 
equations. 

The dynamical degrees of freedom are contained within a and hij. In de Sitter space there 
are no dynamical scalar perturbations of the metric; the "time translation" symmetry implies 
that the usual scalar curvature perturbation C, which we will return to in the inflationary case, 
can be gauged away. (In more general inflationary spacetimes, this can combine with other 
scalars to yield a physical degree of freedom.) In particular, one can choose a gauge where the 
spatial metric is parametrized by 

where 



hij = a\t){e^),, , (2.3) 



det{e') = 1 . (2.4) 

The background equations of motion imply a{t) = aoexp[Ht), where H^ = A/3. 

In an inflating universe, a leading order description of fluctuations is in terms of pairs 
of quanta that are produced at the horizon scale, and then pulled apart by the expansion. 
When we take into account self interactions - either in the matter lagrangian, or due to gravity, 
modes that are produced earlier can influence the correlators of those produced later. A detailed 
calculation of this effect can be done via loop corrections in the in-in formalism (see section 
[3]), but an important point is that the leading IR behavior of this effect can be derived by a 
simple semiclassical procedure. The basic idea, which is illustrated in study of the corresponding 
semiclassical "consistency relations" |9|jll] (and has antecedents in the work of Starobinsky 



12 



and Tsamis and Woodard [13]) is that the early modes are stretched to long wavelength, and the 
later modes therefore view them as a constant background. One can compute the correlators 
of the late modes in the "background" of the early modes, and then do the required average 
over fluctuations of this background to determine the resulting corrections to the correlators. 



To illustrate this procedure (for a review of known semiclassical relations and further detailed 
derivation see appendices A and B) consider the variation of the two-point function, (afciCrfca) 
(where ki are comoving momenta), due to the effect of long wavelength graviton tensor modes. 
This correlator depends on the metric via the magnitude of k. Specifically, for a constant 
background metric perturbation 7jj, we have 



fc — rCiki y Ki\C Jij'^j — rCjfcj 'Jij^ikj + ^'^il'^lj'^iKj + 



(2.5) 



here we use the convention that repeated lower indices are contracted with 5ij. We can then 
expand to quadratic order. 



{(^ki(^k2)^ — ((^fel 0^^2)0 



+ 



d 

-lijkikj + -^iaijkikj ) — {(TkjO-k^) 



[lijkikj 



d 



{(^k^crki) 



(2.6) 



+ ... 



where the zero subscript denotes vanishing 7. The effect of the early, soft gravitons on the 
correlator can then be estimated by averaging this expression over all such graviton modes; this 
can be though of as averaging over a "large box" (here we drop the zero subscript) 



(o-fciO-fe)- 



(o-feiO-fe) 



(2.7) 



d 



+ ijkikjirnaij)— (o-fc^o-fc^) + -kikjkkki{-fij-fki) . ^^2 



d 



(cTfeiO-fca) . 



To compute the average over tensor modes, we first expand 



7*j(a;) = Yl 



d^k 

(2^ 



bleUkh,it) + b\et;i-khm 



,4 



^ikx 



(2.8) 



Here b^ is an annihilation operator, corresponding to helicity s. The polarization tensors e?- are 
chosen to satisfy the transversality and tracelessness conditions e|j(k) = kie^jik) = 0, together 
with a completeness relation obtained by tracing over spatial indices, e*(k)e*J (k) = 26ss'- The 
mode functions are the same as those for a scalar, Uk, up to normalization; they are most easily 
written in terms of the conformal time 



as 



-fkiv) = V2Ukiv) 



7] = -l/Ha{t) 
H 



1 + ikifjc 



ikrj 



These give a two point function 



{l^A^)lkl{x)) = E / (^^(1 + 9'^')4(q)4I(q 



(27r)3 g3 
4 



(2.9) 
(2.10) 

(2.11) 



and the sum over polarizations is 20 



ujij^kiiq) = '^e'j{q)el*i{c[) = 6ik6ji + SuSjk - 5ij5ki (2.12) 

s 

+ SijqkQi + 6kiqiqj - ^tkQjQi - duqjQk - 6jkqiqi - 6jiqiqk + mj^qi 
where q is the unit vector in the direction q. We will also define 



(7^ (a;)) = -Aiij{x)iij{x)) 



d^q H"^ , „ ., f H 

[iTiy q-^ \lTi 



— 1 + g f] 



(2.13) 



The two different contributions entering (2.7) easily follow from these equations: taking 9 
to be the angle between q and k gives 



-j^iliilij) 



H^ 



(2vr) 



'^\l + q'v')ld9{2sin'9) = ^{j\x)) 



Q 



and 






ilijlki) 



H' 



(2vr) 



dq , 



(1 + gV) I d9{sm' 9) = ^{^\x)) . 



(2.14) 



(2.15) 



The momentum integrals are IR divergent. The leading IR dependence can be parameterized 
by introducing an IR cutoff A/r, and for a given time t*, range up to the value q ~ a^:H 
corresponding to modes just exiting the horizon scale at that time. A physical origin for such 



an IR cutoff arises 21 if inflation began at a time t,, in which case q > A/^j = Ojif . (We may 
take modes with q < Am to be in a different state that is less singular than the Bunch-Davies 
vacuum.) Thus these integrals take the formn 



<7^(-)>. 



H^ 



f27r) 



a,H 



dq 



H' 



.H q 



(2vr) 



log 



a* 



(277)2^^* ^'' 



Jj2 



(2vr) 



The result is that the average (2.7) yields 



((^fci^fc.).) = |l + ^(7'(x)>. 



5 \dk^ 



+ k' 



d 
9P 



(o-fciCTfca) • 



(2.16) 



(2.17) 



where t^, is the time of horizon crossing, given by A; = a*if , for the mode in question. The 
massless scalar two point (Wightman) function 



{^kAviykAV2)) = {27rf5\h + h)UkAvi)UUri2) , 
gives a scale invariant spectrum at late times 771 = ri2, 



{ak,ak,)^i2ny5\h + k2) 



2fc3 



(2.18) 



(2.19) 



In eq. (2.161 a^H appears only as a UV cutoff for our semiclassical relations. 



Precisely for this spectrum, the leading IR correction in (2.17) cancelsn However, other corre- 
lators have different power dependence. An example is 



{&k,&k,)^{2nf6''{h + k2)^k 



(2.20) 



where dot denotes t derivative. So, the preceding steps applied to this correlator yield a 

correction 

4 



1 + I^(7^(^)>. 



(2.21) 



Likewise, when we later treat the case of slow-roll, with a spectrum that is not precisely scale 
invariant, there will be similar corrections. Note that these loop corrections to correlators 
become 0{1) on a time scale t ~ 1/H^ ~ RS after the start of inflation, where R and S denote 
the de Sitter radius and entropy, respectively. 

Analogous arguments apply to the graviton two-point function, with the result 



(7fci7fc2)^ 



l + ^(7^(-)>. 



5 \dk^ 



+ k' 



d 
9P 



{ik^lk;) 



(2.22) 



Scale invariance of the graviton correlator in de Sitter likewise implies a vanishing leading 
correction in this case, but this can be altered by slow-roll. 



3 In-in calculation of leading one-loop corrections 

The semiclassical analysis of the preceding section may seem heuristic, and it is important to 
check the method via an exact calculation. What is needed is the one-loop correction to the 
two-point functions in the "large box" using the full Schwinger-Keldysh or in-in formalism. 



As in |9|, this calculation is set up by inserting the the ADM decomposition of (2.2) into 



the action (2.1), which gives the lagrangian 



£ 



a 
~2 



+N-\a - N'd.af - Na-^[e-^Y^di(jdja] 



(3.1) 



where Eij = \{hij — ^iNj — VjA^j) is the rescaled extrinsic curvature and R^^^ is the curvature 
scalar of the three- metric hij. Since the lapse A^ and the shift N^ are lagrange multipliers, they 
were eliminated by their equations of motion, working order-by-order in perturbations about 
dS, in |9]. 

To lowest order we have A^ = 1, Aj = and obtain the free action of gravitons 7, 



S, 



1 



d?x dt a^ [iijiij - a ^dkjijdk'jij] 



(3.2) 



■^This cancellation, which preserves scale invariance, can also be verified to hold for a massless scalar in D 
dimensions, using (cr^) ex k^~^ and the statement {jijjki) oc 2dijdki — [D — l){5ik5.ji + SuSkj)- 



and scalar s a. 



So 



/ d^x dt a^ [d^ - a "^diadia] 



(3.3) 



The free fields derived from these actions satisfy the free equations of motion and can be 



decomposed into the mode functions given in (2.10). With scalar creation and annihilation 
operators a|^, Ok, the the scalar mode decomposition is 



crlx 



d^k 

(2^ 



a^Ukit) + al^U;{t) 



^ikx 



(3.4) 



and the graviton mode decomposition was already given in (2.8). 



cubic order, we have 



We focus on the scalar correlators {o'ki{ilo)'^k2{Vo)) ^^'^ (^fci('7o)^A:2('7o))- The interactions 

re obtainec 

--fijdicrdjo- 



contributing to their one- loop corrections are obtained from the second line of eq.(3.1). At 

a 



C, 



At quartic order, the lagrangian is 

£4 



-nulijdiadja + 



(3.5) 



(3.6) 



where terms with more derivatives are suppressed. These terms arise from the order-by-order 
elimination of the lapse and shift, and are given in |9 11 23]. We find that the leading one- loop 
IR contributions come only from the interactions (3.5), ( 3.6[ )j^ 

In the in-in formalism the expectation value of any operator O (evaluated at time r/o) is 
given by 

(fi| C(r7o) \n) = (0| T ('e^/-°«='^'?^^') 0{r]o)T U-'I-^^r^HA |g^ (37) 

where \Vt) is the vacuum of the interacting theory, |0) is the vacuum of the free theory, T 
and T are time ordering and anti-ordering operators, respectively, and Hi is the interaction 
Hamiltonian for time rj. Two topologies of diagrams give leading IR contributions to scalar 



two point functions at one loop. The first arises from expanding (3.7) to quadratic order in the 
cubic interaction 



H3 = -- d^x a^'jijdiadja 



{3.t 



and is pictured in fig. 1. The second comes from expanding (3.7) to linear order in the analogous 
H4, and is drawn in fig. 2. 

Two kinds of propagator enter these expressions. The first are the Wightman functions. 



The scalar Wightman function was given in (2.18) 



MvWW)) = i2nY6%k + k')Wk{v,v' 



(3.9) 



'''This is also checked in '23j, although that reference obtains different results for the contributions that we 
calculate below. 




Figure 1: Correction to the scalar two point function from an intermediate scalar and graviton. 




Figure 2: Correction to the scalar two point function from a graviton bubble. 



with 

The tensor Wightman function is 



with UJij^kiik) given in (2.12). One also needs the Feynman functions, 



Gfc(r/, r/') = e^'n - n')Wk{v. V') + 0{ri' - r^)Wk{ri\ 77), 



(3.10) 
(3.11) 

(3.12) 



and correspondingly for the tensors; complex conjugates of these expressions also enter the 
diagrammatics. 

There are two types of contributions from the diagrams of fig. 1, and one from the diagram 
of fig. 2. Let us parametrize their contributions to the two point function in terms of quantities 
A, _B, and C respectively, as 



{akijDak^iji)) = {akiv)(^k'iv))o + {2nyS%k + k')[Ak{r]) + Bk{r]) + Ck{v) 



(3.13) 



leaving out contributions of other diagrams with more derivatives, which as we have noted are 
subdominant in the JR. One can straightforwardly find the first and second terms from the 



expansion of (3.7) in H-y. 
Akiv) = -4Re ^ "^'^ 



Bkirj) = 4Re 



(27r) 
(27r)3 



n 



dr]i 



n 



\ {Hiii) 



drji 



j(^ij,kii^)^i^j^kkiGg{r]i, r]2)Gk-q{r]i,V2)Wk{r], ViWkiv, V2) 

(3.14) 

yyijMi'i)^i^j^kkiWq{rji,rj2)Wk-q{r]i,r]2)Wk{r]i,r])Wk{r],r]2) . 

(3.15) 

In deriving these expressions we use transversality oiuij^kiifl), and can also use Uij^ki{'i)kikjkkki = 
k^ sin^ 9, where 9 is the angle between q and k. Likewise, the third term arises from the ex- 
pansion in if 4: 



\iHV^ 



Ckiv) = Re 



-2i) 



d?q 

(2^ 



di 



/\2 



(Hv') 



—uJimjkikjWg{ri',r]'), Wk{v,v') 



(3.16) 



in which we can make the replacement uu^ijkikj = 2k^sm^ 9. If we instead want contributions 
to the two point function of &k, we differentiate both occurrences of r] in these expressions with 



respect to t. From (3.10) and (2.10), this corresponds to the replacement 

{1 + ikr]) -^ -H{krif 



(3.17) 



in each of these expressions. 



Let us begin with Ak, which accounting for the time ordering becomes 



Ak{r]) = -8Re 



dr]i 



dr]2 



d?q 



;k^ sin^ 9Wg{r]i, r]2)Wk-q{r]i,r]2)Wk{r], r]i)Wk{v, V2) 

(3.18) 
In this, we substitute (3.10) and (2.10). The integrals over 772 and rji may then be performed 



exactly, resulting in a somewhat complicated function of the momenta k, q, and p = k — q. 

To do the angular integral, we use 



p^ = /c^ + g^ — 2kq cos 9 



(3.19) 



givmg 



sinede 



pdp 
kq 



(3.20) 



to change variables to an integral over p and q, with p ranging from \k — q\ to k + q. The 
integral over q, which is UV and IR divergent, may be broken up into the two ranges q < k and 
q > k, with an explicit IR cutoff Am used on the lower range, and Auv bounding the upper 
range. (There is no divergence at g = A;.) Then, the resulting expression can be expanded in 
Ajji, 1/AuY, and ki], with small kr] corresponding to large scales/late times. 

We do not show the intermediate steps, as they are algebraically complicated. Analogous 



steps may be performed on the expression (3.15). Combining these results gives 



Akiv) + Bkiv) 



2A?3(27r)2 



2 
3 



2 \og{k/AiR) + ^ log(A,;v'/A;) + ^ + 0(fcV, A?^, l/Aly] 



101 
"90" 



(3.21) 

where one can compute the higher-order terms explicitly (see appendix C). The third integral 
Ck is more straightforwardly performed, and has expansion 



CM 



2 Jt2 



2A?3(27r)2 



[2 log(A,^,/A^v) + 0{kW, AJR, 1/A 



2 

UV 



(3.22) 



Combining these, we find that the IR divergent corrections cancel, as we found in ( 2.17[ )j^ 
Likewise, repeating the steps for {ak^ak2), using the replacement (3.17) gives 



(o-fciCTfca) ~ (o^fciCrfc2)o 



^-]|(^'°«"''«/'=' 



(3.23) 



which, using k = a^,H and (2.16), is in perfect agreement with (2.21) 



Thus, the results of this section are a non-trivial check of the validity of the semiclassical 
relations of the preceding section for finding the leading IR corrections to correlators. 



^ After the first version of this paper appeared, it was shown in 



22 



that in the full expression, taking into 
account other UV divergent terms, the log k dependence cancels, leaving log(A/;^;/At/y). Comparing with (2.16 1, 
a^H is thus replaced by the true UV cutoff h^uv- 



10 



4 Slow-roll 



Having performed non-trivial checks on the vahdity of the semiclassical relations in de Sitter, we 
now extend their use to slow-roll inflation. We assume that the methods given in the previous 
section can be straightforwardly extended to the slow-roll case, and perform various consistency 
checks on this assumption. For definiteness and simplicity we study the example of single field 
slow-roll inflation, although most of the subsequent discussion can be trivially generalized to 
other models of inflation. 

In single field slow-roll inflation, inflation is driven by a scalar field with a very flat but 
otherwise arbitrary potential V{(f)), 



S = lJv^[R-d, 



As usual, we define the slow-roll parameters 



1 rv^ 

V 



V 



- 2V{<P)] 






(4.1) 



(4.2) 



Perturbations are again conveniently analyzed using the ADM decomposition (2.2). Now 



however, the scalar degree of freedom mixes non-trivially with the previously gauge scalar degree 
of freedom of the metric. Different gauges may be used, but for the following considerations 
we find it useful to work in the comoving gauge, where the fluctuations in the inflation vanish. 



= 0o(t) , 

with 00 the slow-roll solution, and the spatial metric is parameterized as 



hi 



a\t)e'He^)., . 



(4.3) 



(4.4) 



Here again det{e^) = 1, and in this gauge ( can be non-perturbatively identified with the gauge 
invariant curvature perturbation which is conserved on large scales (9,24,25 



Now, both scalar and tensor fluctuations couple, and contribute to loop corrections. We 
consider them in turn. 



4.1 Tensor loops 

The analysis of tensor loop corrections to scalar propagators parallels closely that of section [2j 
Specifically, we first begin with 



(CfciCte)' 



(CfclCfc2)f 



+ 



''JijI^iKj + 'Jil'Jljkikj 



d 



iCkiCkz) 



+ 2 (^«i^«%) 



d 



iCkiCki) 



(4.5) 



+ ... 
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We then average this over tensor modes in the "large box," following the steps of section [2| 
with the result 



(CfciCfc2)' 



l + 3(7^(^)>. 



5 \dk^ 



+ k' 



dk^ 



(CfciCfca) 



(4.6) 



Now, we have the IR divergent integral 



f(^)>. 



a^H 



dq H^ 



(4.7) 



with time-dependent (q-dependent) H. Without scale invariance, the correction doesn't cancel. 



If Ug is the scalar spectral index (n^ = 1 as usual for scale invariance), the expression (4.6) 
becomes 



(CfciCfc2)^ 



Likewise, for tensors we find 

(7fci7fc2)- 



5 






(CfciCfca) 



(4.8) 



(7fci7fc2) 



(4.9) 



where rit is the tensor index. When applying {d/dk'^Y, we will in general also get a contribution 
from the running of the spectral index dng^t/dk'^. Above we have however neglected the scale 
dependence of the spectral indices to leading order in slow-roll. More detailed expressions are 
given in appendix B. 



4.2 Scalar loops 

Accounting for scalar loops follows a similar procedure. A long-wavelength background scalar 
( produces the modification 

k"" ^kl = {e-~^kf . (4.10) 

Then, taking into account the dependence of the momentum-space measure on (, we have 



(CfclCfc2)i 






''^"(C(e-%)C(e-%)) 



(4.11) 



(CfciCfc2)o- (^s-1) C(CfciCfc2)|o+ ini^s- l)^ + «s) CC(CfciCfc2)|o + 



Here we have introduced the running of the spectral index j^a^ = dns/dln{k), and we only give 
expressions to leading order in slow-roll. More complete expressions are given in appendix B. 

^In the first eprint version of this paper, we had neglected the running of the spectral indices, but it was 



subsequently pointed out in 26 that these are generically relevant in the case of scalar loops. 
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So, then again we obtain the effect of the soft scalar mode on the two point function, by 
averaging over the two-point correlation function in the large box 

{{CkA.)c) ^ (Ck-Xk,)o + (^K - 1)' + as) (a,a.)o {c'{x))^ , (4.12) 

with 

Similarly, introducing at = dnt/dln{k), 

{{ik.lk,)^) ^ {lk,-fk,)o + {^{ntY + at) {ik,lk,)o {C\x))^ , (4.14) 

here we use the standard convention where for historical reasons n^ = 1 corresponds to a scale 
invariant scalar spectrum, while Ut = corresponds to a scale invariant tensor spectrum. The 
corrections are suppressed by slow-roll factors 

ns-l=2r]-Qe , rit = -2e ; (4.15) 

they vanish in pure de Sitter, where the scalar mode ( decouples, as expected. 

4.3 Alternative explicit check 

As a check on the semiclassical relation with scalars, we can use an explicit calculation of the 
second order gauge invariant curvature perturbation |27|j7^''^''. Ref. 27 obtained in the case 
of as = 

7^(') = -^(r^,-l)(7^«)2+X (4.16) 

where the non-local contribution X is determined by the integral 



-- hi^^^^fdt + (e - 77)A-la,7^(l)9*7^(^) , (4.17) 

and on large scales z/'^^^ = tRS^\ The one-loop contribution to {IZkR-k') is then given by 



TZfTZ^^) 



First we note that on large scales all higher order gradient terms and time derivatives of 
if)^^^ vanishes to leading order in slow-roll. So on large scales to leading order, we have 

X ~ (e - r/)A-lai7^W9*7^W, (4.18) 



^The curvature perturbation TZ is related to our C, hy TZ — ~Q &s discussed in the second footnote of 29] 



around eq.(3). The curvature perturbation in 27 is related to the conserved curvature perturbation TZ only in 



the case of a^ — 0, as was pointed out in 28 . We thank Filippo Vernizzi for bringing this to our attention. 
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which in Fourier space gives 

IW~(^-'rt/(gs'^4^K<')<, (4.19) 

SO the potential IR divergence in 7?.q from the lower limit of the g-integral, is always regulated 
by a multiplicative power of q in the limit g — )■ 0, and there are no IR divergences from this 
piece to leading order in slow-roll. 

So, let us finally look at the contribution to l/Tci. 1^,1 ) from the local piece. From that we 

obtain 



- ^(--l)V||/||^'^«>('^L\<-«). (4.20) 



If we split the integral over q into an IR part Ojifj < q < k = a^:H^: and a UV part q > a^H^, 
then from the IR part of the integral, we obtain to leading order in slow-roll 

Kfwg')" = l(n, - I? (k™<) (k'"'(-)). , (4.21) 



in agreement with eq. (4.12) in the case a^ = 0. 



4.4 Comments on relation to SN approach 

In the 6N approach, the relative curvature perturbation on super-horizon scales can be de- 
scribed by the evolution of causally disconnected regions, which are assumed to evolve like 
separate unperturbed universes. In this approach, which corresponds to the leading order in 
the gradient expansion, one has ( = 6N, where A^ is the number of e-folds of expansion between 



an initially flat spatial hypersurface and a final uniform density surface 25 , 30 33 . This im- 



plies that in single field slow-roll inflation, one can write the comoving curvature perturbation 



as an expansion of the field perturbation at horizon crossing 34 



C = A^'¥Uo..e.o.. + I N"6<p ^ <50|,„,,_. + . . . (4.22) 

which now describes the classical relation on super-horizon scales between ( and 50, or equiva- 
lently it can be viewed as the super-horizon part of the gauge transformation between the two 



variables in uniform curvature gauge and comoving gauge respectively. Note that in eq.(4.22) 
above, prime denotes derivation with respect to the field 0, and star denotes a convolution of 
momenta in Fourier space. As an example, by use of the chain rule, one can to linear order 



'For more details, see in particular the discussion in section IV of 25 
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write A^' = dln{a)/d(f) = {dt/d(f)){dln{a)/dt) = H/(p, in order to recover the linear relation 
between ( and 50 on large scales. 

In this approach one can calculate the one-loop corrections to correlation functions of (, 
by calculating one-loop correlation functions of 6(j) at horizon crossing, and relating them by 
the 6N approach, which is equivalent to the appropriate gauge transformation. In Fourier 



space, the 6N transformation in eq.(4.22) will at the non-linear level introduce new momentum 
integrals, called 6N loops in the literature, which can contain additional IR divergences. 

However, it is known that already at horizon crossing the correlation functions of 6(f) contains 



IR divergences 35 36 , which cannot be captured by 6N loops. Thus the IR effects computed 



from 5N loops I34l37 40 are incomplete. One important manifestation of this incompleteness 
is that the tensor modes do not transform in the 6N approach (they are gauge invariant on 
super-horizon scales) and hence there are no SN loop corrections for correlation functions of 
tensor modes, while the calculations of the previous sections clearly show that the correlation 
functions of tensor modes receive IR one-loop corrections. 

5 Contributions to cosmological observables 

The leading order IR loop corrections that we have found can produce corrections to cosmolog- 
ical observables, which can become significant for sufficiently long periods of inflation. In this 
section we first discuss corrections to the tensor/scalar ratio and nongaussianity parameters, 
averaged over large distances, and then estimate the size of the corrections in a model where 
slow-roll begins with an exit from the chaotic regime. 

5.1 Tensor-scalar ratio 

The tensor-scalar ratio is one of the main observables in infiationary cosmology. The absence of 
a detection of primordial tensor modes severely constrains the tensor-scalar relation and leads 
to an exclusion of A0^ infiation (see fig. (14) of ^IJ). In slow-roll infiation there is at tree-level a 



definite prediction for the tensor-scalar ratio, called the tensor-scalar consistency relation 42 
We show that the IR loop corrections that we have found shift this relation in the "large box" . 
The tensor-scalar ratio, r, is defined by 

r = 4i^ (51) 

where at tree-level in slow-roll we have (|CfcP)o — (l/^e) (l7fcP)o' which implies the single field 
tensor-scalar consistency relation (see e.g. eq.(23) of |42| ) 

r = 16e . (5.2) 



Now, let us consider the one-loop result. For the scalar spectrum, we combine (4.12) and 
( |48l ) to find 



((CfciCfea)) - (CfciCfe2)c 



1 + (^K - 1)^ + «.) (C^(x)>, + \^\^ (7^(-)>. 
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(5.3) 



Likewise, for the tensor spectrum, we combine (4.14) and (4.9): 
((7fci7fc2)) = (7fci7fc2)o 



i + (^M^ + «0(C^(-)>. + \^|(7^(-)>. 



(5.4) 



Using these gives the corrections to the ratio at one-loop order: 



16e 



1 + 



^2 + 2at - {{us - If + 2a,) 



(C^(^))* + 



{rit - ?>)nt - {us - 4)(n, - 1) 
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(f(^))* 



(5.5) 

From this we see that the loop corrections we have found generically alter the slow-roll tensor- 
scalar consistency relation. So while the loops don't change the scaling of the two-point func- 
tions, they do contribute shifts to the averaged tensor/scalar ratio, as anticipated in [36]. In 



cases where the scalar fluctuations dominate (for one such case, see section 5.3), note that 
nf + 2at > {us — 1)2 + 2as would imply r > 16e. 



5.2 Non-Gaussianity 

The non-gaussianity of primordial cosmological perturbations arising from the three point func- 
tion of curvature perturbations is usually parametrized by the dimensionless parameter /tvl, 
which gives the relative strength between the three point function and the two point function. 
Just as the power spectrum, P({k), of curvature perturbations is defined by 



(Ck.Ck2) = (27r)^5^(ki + k2)Pc(^) 



(5.6) 



which is related to the usual scale independent power spectrum by P({k) = {27i^)V(;{k)/k^, one 
can similarly one can define the bispectrum 5^ by 



(C(ki)C(k2)C(k3)) = i2nf5iY,^a)Bdk,MM) , 



(5.7) 



The strength of the non-Gaussian signal in the bispectrum can conveniently parameterised in 
terms of the dimensionless non-linearity parameter Jnl, by defining 



6 

5^ = —-fNL[P(iki)P(^{k2) + 2 permutations] 





(5.^ 



Now, from the semi classical relations, we can also compute the effect of the long wavelength 
background of scalar modes on the three point functions 



{CkiCk2Cks)c 



. 7 9 1-2 -9' 



-9f(C(e-^"A;0C(e-%)C(e-%)) 



(5.9) 



(CfciCfc2Cfc3)o - (^3 - 1) C (CfeiCfc2CA:3)|o "^ 7) (CfciCfc2Cfc3)lo + • • • 
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where like in 43 , we define 



ns - 1 = 2(n, - 1) + 



dlog{f, 



NL 



d\og{k) 



(5.10) 



Then again completely analogous to the calculation of the effects on the two-point function, 
we obtain the effect of the soft scalar mode on the three point function, by averaging over the 
three-point correlation function in the large box 



iCkiCkiCks), 



{(kA.(ks)o + lin,- 1)2 (^,(^^3)0 (C'(^)>, • (5.11) 



For simplicity, let us consider the squeezed limit ki << k2, k^ of single field slow-roll 
inflation, where at tree-level /at^ = (5/12) (rig — 1) |9|. In this case we then have 



Br 



[n. 



1) 



l + 2(^3-l)^(C'(x)>J 



[P^P{ki)PP{k2) + 2 permutations] 



(5.12) 



where P^\k) is the tree-level power-spectrum. By using the formula for the two point function 



in eq. (4.12) 



P^{k) = P^'\k) 



i+(^K-i)' + «.)(C'(^)>. 



(5.13) 



we obtain in the squeezed limit 

5 



/ 



NL 



12 
5 

12 



-^{n^-lf-{{n,-lf + 2as) 



(n, - 1) 1 1 

(n,-l)[l + ((n,-l)2-2a,)(C2(a;)>J , 



c^(-)>. 



(5.14) 



where we in the last line used that in the squeezed limit Jnl is independent of k. This indicates 
that at one-loop IR effects shift the magnitude of the non-Gaussianity averaged over the "large 
box" upward in single-field slow-roll inflation. We expect tensor modes to also contribute cor- 
rections, although in the examples considered in the next section these are relatively suppressed. 
We also expect corrections for more general momenta, and higher-order correlators. 

5.3 Example: chaotic inflation 

In chaotic inflation the initial conditions for the inflaton field are determined by large quantum 
fluctuations in the "self-reproduction" regime. In this regime the quantum fluctuations of the 
inflaton field are larger than the classical evolution of the background field, 50 > 0cAt (with 
At ~ ^/H), and the quantum fluctuations randomly cause some regions to inflate, while in 
other regions the field value set by the random quantum fluctuations does not allow for inflation 
to begin. 
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A simple model is provided by a monomial inflaton potential 

V{<P) = XM'f^Y ; (5.15) 



•■p 



inflation takes place when e < 1, which implies large field values (J) > (pf = a/\/2. At very large 
field values, when the potential energy density becomes of order 0{Mp), the low energy effective 
field theory approach to gravity breaks down. This happens a field values > 0p = \~^^°'. 
However, the regular perturbative description of infiaton quantum fiuctuations on a classical 
background evolution only becomes unproblematic at a lower intermediate energy scale, given 
by the self-reproduction scale, at field values 

< 0, = (127r'aVA)^/("+2) . (5.16) 

Thus, the conventional infiationary regime is determined by 0j > > 0/. 

In chaotic infiation, the total infiated volume can typically be much larger than the ob- 
served universe, if infiation is assumed to begin just at the end of the self-reproduction regime. 



Therefore, in such a scenario we might expect large IR effects. We see from eqs. (5.3), (5.4) 



and (5.12) that the IR effects of the scalar and tensor loops provide fractional corrections given 



by the quantities 

9IR = {\{ns - If + as) {C\x))^ , hjn = (n, - 1) {i\x))^ , (5.17) 

evaluated at the time t* where the visible modes exit the horizon. Thus, in order to evaluate 



their size, we need to evaluate these quantities for the slow-roll potential (5.15). 



The scalar variance is given in eq. (4.13), and can be computed in terms of the change in 
the field value and the parameters of the potential. The build up of the classical perturbation 
({x) on large scales can be regarded as Brownian motion of ({x). The classical perturbation 
C(x) receives a kick of size {l/y2e)H/27i in a random direction every time interval At = 1/H, 
from the quantum mode (q with q = aH, which exits the horizon and becomes classical in the 
given time interval. In de Sitter, with constant H, the mean squared distance to the origin 
is proportional to the time elapsed [8], i.e. (o"^(x)) = H^t/^ATx"^). In slow-roll, since H are 
larger at the beginning, the kicks are larger at the beginning, and the variance is expected be 
larger than H^ {t)t / {Att'^) . In order to compute the variance of ( during infiation, it is therefore 
convenient to express the integral over q as an integral over field values at time q = aH, when 
the given mode is exiting the horizon. 



Specifically, to leading order in slow-roll 

d log q d log q d log a dt H 3H^ 



d(f) d\oga dt d(j) (j) V a 

Then, using H^ = V/3 and (4.2) gives 



(5.1^ 
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Note that this is dominated at the upper value, ~ 0i, so we will neglect 0^,. Next, the 



condition for the boundary of self-reproduction (5.16) gives the value of (pi in terms of a and 
A, and thus 

{eix)l ^ ^^^{l2.W/Xr^--^^ , (5.20) 



The coupling constant A is then fixed by requiring that the observed value from WMAP 7 [44] , 
P^(fc) = (2.43 ± 0.11) X 10"^ is obtained. This gives 

A = 127r2a2(2aAr,)-i-"/2p^(A;) (5.21) 

where A^* is the number of efolds left of inflation when the visible modes exit the horizon. 



N, = j Hdtc^J _rf0=— (0^-0^). (5.22) 

Depending slightly on the reheat temperature, we can assume that when the observable modes 
exit the horizon N^ = 60 in order for inflation to have time to solve the horizon and flatness 
problems of the standard Big Bang model. Thus, we rewrite 

<C^(x)>^ - ^n(^)-^/^"^^^ • (5.23) 

We can also express the spectral index n^ — 1 = 2?7 — 6e in terms of A^*, using 

^ ^ (5.24) 



1\V ) 1(^1 4Ar, 



These give 

Similarly one can compute 



V = — = ^ — - = • (5.25) 



ns-l = -{l + a/2)/N, . (5.26) 



as = -{l + a/2)/N^ . (5.27) 
Combining these results, we find the fractional size of the scalar contributions, 

m = (^(n. - D' + a.) (f (.». :. {^^^^ (n(»:))-""-''' ^ (5.28) 

Taking N^ = 60 and V^{k) = (2.43 ± 0.11) x 10"^ [44], the fractional size of the scalar 
correction parameter is to first order in slow-roll 

^ gjR ~ -1383 (5.29) 

^ giR^O (5.30) 

^ 9iR^8 (5.31) 

^ 9iR^5 (5.32) 

19 



a = 


= 1 


a = 


= 2 


a = 


= 3 


a = 


= 4 



We see that at first order in slow roll, there is a cancellation for the case a = 2, and the leading 
contribution conies from higher-order in this expansion. 

The tensor correction parameter is calculated similarly. We find 

<7^(^)>. = r . \,, (€^' - €^') > (5-33) 

and so 

hjn = in,-l){^'ix))^ ^ — . (5.34) 

The tensor contribution is thus suppressed relative to that of the scalar. For general cases 
one finds gjn > 1 and the corresponding perturbative corrections become large before inflation 
ends, if it is assumed to begin just below the self-reproduction regime. Thus, correspondingly, 
the corrections to the tensor/scalar ratio and nongaussianity parameters, averaged over the 
"large box," are large. 

6 Comments on IR effects 

Infrared divergences in the propagator and variance of a light field have long been known and 
explored in the inflationary literature. We have found that for long periods of inflation, or 
pure de Sitter space, these divergences are associated with large loop corrections to physical 
quantities in the total inflated volume ("large box"). These can be thought of as due to 
accumulation of long wavelength fluctuations, that provide significant corrections to correlators, 
and this phenomenon can be directly analyzed in a semiclassical approach, using so-called 
"consistency relations." 

One might at first sight be puzzled by the appearance of large logarithms like we have found. 



given the analysis of 46 , which appears to argue that large logarithms don't arise in correct 



calculations, and in particular are forbidden by the symmetry a — )■ Xa, x — ?■ x/X. However, the 



analysis of 46 applies to logarithms that are UV large, and not to IR-divergent logs. Indeed, 



note that the IR-regulated expressions, e.g. (2.16), respect this symmetry. 



In a future publication |47j we will explore other manifestations of these and related IR 
effects. Such IR effects are clearly physical. For example, accumulation of IR scalar fluctuations 
serves as the driving mechanism for self-reproduction. Here we have shown that a similar 
buildup of graviton fluctuations leads to corrections to correlators that may be significant. 
Such strong IR behavior also enters other physical quantities, such as the tri-spectrum. For 
such IR divergences in modes coupled with gravitational strength, the characteristic time scale 
at which the corrections become strong is t ~ RS, where R and S are the de Sitter radius and 
entropy, respectively. 

The question of how to handle such IR divergences has received considerable discussion. For 
example, one might try to resum the IR effects to give a corrected background, and then expand 



about this background. As a first approach, consider for simplicity eq. (4.11) and eq. (4.12) 



although the argument will be exactly the same for the tensor modes in eq. (4.6) apart from 
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the trivial extra complication of polarization sums. If we continued the Taylor expasion in 

d 



eq. (4.11) to arbitrary order, eq. (4.12) would become 



(CfciCfca) 



(a.a.)o + (2vrf<^'(ki + k2)J] ^2^^, 



^^^-^^"(c^(.)):(3 + 



d\og{k) 



2n 



{CkiCk2)i 



CXJ y 

ra=l ^ 



ris-l 



+ a/2 



2n 



iCix)): 



Pdk) 



(6.i: 



where L'^{x) is the generalized Laguerre polynomial, and in the second step we restricted 



ourselves to monomial potentials of the type eq. (5.15 ). Interestingly this sum does not converge 



in general, and the perturbative treatment seems to break down completely as the variance 
becomes larger at long distance, and all terms contribute to the sum. Since {5(fP'{x)) ~ H^t, 
we see this happen past t ~ RS. For example, in the model with a = 2 studied in section 5.3 
classical evolution led to {ug — lY {C,"^ [x)) ^^ = 450. If we were to take the resummation at face 

value, it would give a divergent result for / {CkiCk2)c)^ ^^ ^^^ ^^^ '^^ chaotic inflation, indicating 
that we had entered a non-perturbative regimq^ with quantum fluctuations overwhelming the 
classical background evolution in the "large box" . While elsewhere it has been advocated W^ 
that loop effects destabilize de Sitter space by negating the cosmological constant, the physics 
we apparently encounter here suggests another type of instability of inflationary universes, to 
growth of large fluctuations of fields and geometry on large scales. This is similar to self- 
reproduction, but apparently more general. 

While various approaches to handling this regime have been advocated, for example re- 



summation via a dynamical renormalization group approach 48[ 49 or via stochastic ap- 
proaches 



12 13,50 



no method has yet been found for a precise and complete treatment 
of the quantum dynamics into the regime where such effects become important. It is not clear 
how the large corrections can be absorbed into a background correction. It may be that the 
regime where such effects become important represents a true boundary in our description of 
physics via perturbative gravity, and that complete description requires fully non-perturbative 
gravity, quite possibly with new dynamics. 

While some IR divergences are dimension dependent, note that the ones discussed here ap- 
pear to persist in higher dimensional gravity. In D spacetime dimensions, dimensional analysis 
tells us that in a time 1/H, the fluctuation of a scalar field is 50^ ~ H^~^. This in turn 
indicates growth {6(f)'^{x)) ~ H^^^t, again becoming large on a time scale t ~ RS. 

Note also that IR divergences are typical indicators that one has not asked a sensible 
question, as in the story of soft photon effects in QED and the need to calculate inclusive 
sums. In the regime where IR effects become important in gravity, it may be that there is 
a good local description of the dynamics, as in pictures of one inflating bubble of a chaotic 
universe, described for sufficiently short times. That is, there may be good local questions, 
such as the evolution of quantum fields in a limited region of the chaotic universe for a limited 



^Note that this behavior is different from the behavior of an 0{N) invariant test scalar field in the large N 
limit, where the exact resummed correlation functions approach a constant for large values of gifj 50 . 
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time, and resummation and/or renormalization group methods may be useful for describing 
such a piece of the universe, where for example the fluctuating backgrounds {e.g. the scalar 
vev, and thus vacuum energy in a chaotic scenario) take particular valuesPj For example, we 
might view the correlation functions in the large box as the average correlation function in a 
"small box" (Hubble sized region) within the large box. This interpretation is manifest the 
way we calculate our semiclassical relations. The average correlation function in the large box 
is not directly measurable for a local observer, and some IR effects might be absorbed into the 
uncertainty of the state of the fields of the small box in the large-box ensemble; in particular we 
might not live in an average small box within the large box. This interpretation also resembles 



some of the ideas in 38 , 39 . In particular, in the special case of single field inflation the end 
of inflation is determined by the v.e.v. of the inflaton, and one might condition the probability 
distribution of the correlation functions on only 60 e-folds being left of inflation in the small 
box, using the inflaton as a clock [5l] (see however also [52]). But more global questions, such 
as that of determining the quantum state of the universe on a global slice, and/or at long times, 
and associated quantities such as averages over the slice, may not have answers - at least not 
in perturbative gravity. 

Indeed, this viewpoint that such phenomena could represent a true breakdown requiring 



new nonperturbative dynamics and mechanisms has been previously suggested in 15,16 . In 
fact, there is a situation in black hole physics that is quite similar, and is associated with the 
so-called information paradox. As originally argued by Hawking (l4J, perturbative methods 
imply that black hole evaporation destroys information, but this is very problematic p] A more 
careful argument requires calculation of the quantum state of fields on a spatial slice spanning 
both the interior and exterior of the black hole. Such "nice slices" are analogous to the spatial 
slices of de Sitter space or slow-roll |15] - [l8|j^ Note that while in cosmology one may view the 
state on such a slice as unmeasurable, in the black hole context it is more directly connected 
to measurable quantities since ultimately we imagine measuring the complete quantum state 



of the system {i.e. all outgoing Hawking radiation). Ref. 15,16 argued that a perturbative 
calculation of the state on such a slice breaks down due to accumulated effects of fluctuations, 
on the time scale t ~ RS, where now R and 5* are the black hole radius and entropy]^ We have 
demonstrated an analogous growth of perturbative corrections in cosmology. In the black hole 
context, 15 16 proposed this as a resolution of the information paradox, since such a breakdown 
means that there is no sharp calculation of the state on the slice, and thus no sharp calculation 
of the missing information. While this would eliminate the contradiction producing a paradox, 
it leaves an information problem, which is to determine the non-perturbative mechanics that 
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Moreover, a more basic description of such local questions may be through relational observables, such as 



the proto-local observables of 53 . 

iipor reviews of the information paradox see 



54 



55 



^^Analogies between black holes and de Sitter cosmology have long been pursued; for other examples see 
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e.g. 56 , \IS\ 



Rei. [TTI 18 suggested a breakdown of a different origin due to small spacing between nice slices. Note, 



though, that they also argue that any attempt to regulate inflation through slow-roll will fail after a time scale 
t ~ RS as a transition to eternal inflation will result, which does fit with our discussion. Also, [57 argued that 
fluctuations of observers/observables become important on this time scale. 
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takes over in this regime, and by which the information escapes the black hole, while it is still 
of a macroscopic size comparable to its original radius. We see close parallels between this 
discussion and that of inflationary cosmology. 
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A Review of consistency relations 

In this appendix we will review the consistency relations for the bispectrum of non-Gaussianity 
in the squeezed limit |9, 10 and for the trispectrum of non-Gaussianity from the exchange of a 



long wavelength graviton between two pairs of scalar modes in the counter-coUinear limit 11 . 



A.l Squeezed limit consistency relation for the 3-point function 

Consider the three point correlation function of curvature perturbations (CfciCfc2Cfc3) in the limit 
where one of the momenta is much smaller than the other two ki << ^2, ^s- Then the long 
wavelength mode will rescale the spatial background of the other two modes, as can be seen 
from the form of the metric 

ds^ = -dt'' + e^'^'a\t)dx\ (A.l) 

Now, we can Taylor expand the two-point function of the two short wavelength modes (k2 5 Cfca 
on the background of the long wavelength mode Ci 



d 



+ . . . (A.2) 







and then by correlation of the variation of the two point function of the shifted background 
with the long wave mode, one obtains the three point correlation function in the squeezed limit 

lim (CfciCfc2Cfc3) = (Cfci (CfcaCfca)^ ) 

{ns-l){(k,C-k,){Ck,(ks) ■ (A.3) 



Using the definition of /nl in eq. (5.8), we obtain in this limit from single field slow-roll 
inflation [o] 

Inl = ^(^s - 1) (A.4) 

Note, we could just as well have considered the shift in the background due to a long 
wavelength graviton mode 7^, which would have the effect dx"^ — )■ dx"^ + 'jjjdx^dx^ or in Fourier 
space k"^ -^ k"^ — •yf^kikj. So the we would obtain 
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(CfcaCfca) 



71 



(Cfc2Cfc3)o+ 7li 



d 



dlu 



(CfcaCfca) 



d 



+ ... 



\>>k2>>kzlQ ^'^i^U * i Qjr,2 V^'^aSifca/ 



+ 



where the e^j is the polarization of the Fourier mode 7^^ . One then obtains 

lim (TfciCfcaCfca) = (7fci (Cfe2Cfc3)7i 

In the same way one could calculate 






lim (7fci7fc2 7fc3) 



7fci (7^27^3) 



■3/71 



<J ''^i ^ijl^il^j 



(7fci7-fci) (7^27^3) 



2 fc2 
All of the above relations have been verified by explicit in-in calculations in [9 



(A.5) 



(A.6) 



(A.7) 



A. 2 Counter-collinear limit of graviton exchange diagram 

Consider the contribution to the four point correlation function (CfciCfc2Cfc3Cfc4) from the exchange 
of a graviton between two pairs of scalar modes in the limit of the momentum of the exchanged 
graviton going to zero. Defining ki2 = ^1 + ^2, this limit corresponds to the so called counter- 
collinear limit with ki2 << fci ~ ^2, ^3 ~ ^4- 

In the limit where the momentum, ki2, of the exchanged graviton goes to zero, its effect is 
again to rescale the background dx'^ — )■ dx^ + •yj^dx^dx^ or k'^ ^ k"^ — 'j:^kikj. Taylor expanding 
again on the background of the long wavelength graviton 



iCkiCks). 



d 



+ 



{Ck2Cks)o + 7ii^— B iCkiCks) 

d 

{Ck2(ks)o + IkB^ij^i^i-Q^ (Cfc2Cfe3) 



+ 



{k.i 



one can obtain the semiclassical contribution to the four point function from the correlation 



between a pair of two point functions due to the long wavelength graviton 11 



lim (CfciCfc2Cfc3Cfc4) - {{CkAki) ^. iCkzCki) ^. 
fcl2— >0 \ "°12 "=1 

_ (n^ - 3)2 eijkikj eikhk^ 



(CfciCfca) (7fci2 7-fci2) (Cfc3Cfc4) (A. 9) 



4 A;2 J.2 

This relation has been checked by taking the appropriate limits of the full in-in calculation 



in 11 , where the straightforward generalization to scalar exchange diagrams was also discussed. 
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B Detailed derivations of the semiclassical relations 

In this appendix, we give the detailed derivations of the semiclassical relations for the one-loop 
IR corrections to the two point functions from virtual tensors and scalars in the loop. 

B.l Details of the derivation for virtual gravitons 

We want to calculate (CfciCfc2) o^i the shifted background of infrared gravitons. In the perturbed 
spatial metric 

h,, = a^e^^ [e% ^ a'^'- \h, + 7., + ■ ■ ■ ] , (B.l) 

the effect of an infrared graviton ^^ is therefore to take dx^ — )■ dx^ = dx^ + '-/l^dx'^dx^ + . . . . 
Thus, we want to calculate the average correlation function in the large box, given by 

{{CkA.)B) ■ (B.2) 

This is the average in the large box of {CkiCk2)B^ where (CfciCfc2)B ^^ ^^^ correlation function 
of the short wavelength modes in the background of the long wavelength modes somewhere 
inside the large box. In real space we can then Taylor expand the correlation function on the 
unperturbed background 

(C(a^i)C(a^2))B = (C(a^i)C(a^2))o + 7»fK)^ (C(5i)C(52))|^^=o 

+ ^7.?(a:o)7.1(^o)^^ (C(£i)C(£2))|,,.o + • • • • (B-S) 



Using that 



(^E4(%|W^*'' (B-4) 



we can then compute (CfeiCfe2)B) by a Fourier transformation 

{CkA.)B = J J d'x^d'x2e-'^^-'^^e-''^-'^^ax,)ax,))^ 

= (CfclCfc2)o 

+ [ /"rf33,^^3^^g-m.fcig-iX2-fe2 



(9 1 d d 



d^qi d^q2 




(27r)3 (27r)3 



e^^-'^-e^^-'^- (^-,^2) 



B 



(B.5) 



7s =0 



Since the long wavelength mode 7^ is almost constant over the scales of variation of the 
short wavelength modes, we can choose Xq freely between xi, X2- For simplicity we can take 

^0 = (^1 + :?2)/2. 

25 



In order to evaluate the expression above, it is useful to write 



x^ 



x^x^ 



X^ + -ffjX'x^ + ... 



(B.6) 



and using that momentum transforms inversely, define 



q 



QiQj = q^ - ifjQiQj + ■■■ , 



(B.7) 



such that X ■ q = X ■ q, and 



d dq d 



dlfj dlfj dq^ 



d 



(B.8) 



Note also det[exp(7^)] = 1, so the Jacobian of the transformation g — )• g is trivial, and 5^{qi + 
92) = <^^(^i + ^2)) so we can write 



\SgiS,(?2/ 



B 



{27if 5\q^ + q2) 



1 + ( -ifjm, + ^luiSm + • • • j ^ 



+UiiSmy + ---)^^ + --- 



■^ f)n-^ 



dq^ dq 



{(,(-,) 



(B.9) 



Inserting this into eq. (B.5) and changing the integrations over q into integrations over q, we 



then obtain after taking the average in the large box, and after integrating out three delta 
functions 



((CfciCfca)^) - (CfciCfe2)o 



(B.IO) 



+ 



l{2^mh+h) I ^,j:h.<j 



d 



d d 



4iikBy*{kB)ijhk,j— + e'{kB)ijkikje'*{kB)kihki-^-^ 



Pdk) 



Inserting the expression (2.12) for the polarization sums 



J24ikB)e'*{kB)i,k,k^ = 2k'sm\e) , ^ e^(M,,A;,A;,e^*(MMfcfcA;; = A;Sin^(^) , (B.li: 



and then doing the angular integrals, we finally obtain 



{{CkA.)B) = (^,^2)0 + i2rr)'6'{h + h)- {i\x))^ 



5 \dk^ 



+ k' 



_d_ 
dk^ 



Pdh), (B.12) 



in agreement with eq. (4.6) 
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B.2 Details of derivation for virtual scalars 



In this subsection we want calculate {{CkiCk2)B) from the shifted background of infrared scalar 
modes. For convenience let us neglect for a moment the tensor modes and expand the perturbed 
spatial metric as 

h,j = a'e^^d^j ^ aHi,(l + 2C + . . . ) (B.13) 

The effect of an infrared mode C,^ is therefore to take x — )■ x = x + C,^x + . . . . 

Like in the previous subsection, we Taylor expand the correlation function on the unper- 
turbed background 



{C{X,)C{X2)) 



B 



d ^ ^ 

(C(xi)C(x2))o + Cb{xq)-q— (C(a;i)C(a^2))|^^=o 



1 d d ^ ^ 



Cs=0 



Using that 



C{x) 



(2^ 



Ck{t)e 



ik-x 



we can then compute (CfeiCfe2)B) by a Fourier transformation 




rf^Xirf^X2e-*^i-''^e-*^2-'=2(C(xi)C(x2))B 



(CfciCfc2)o + 
_d d_ 




Sxid^X2e 



-ixi-ki —ix2-k2 




d^ke d^qB^i^,, 




d^qi d^q2 



(27r)3 (27r)3 



(27r)3 (27r)3 
6(;Bjxi-qijx2-q2 /((e~'^^gi)C(e"^^ ^2) ) 



+ ... (B.14) 

(B.15) 

(B.16) 



Cs=0 



Since the long wavelength mode 75 is almost constant over the scales of variation of the 
short wavelength modes, we can choose xq freely between xi, X2- For simplicity we can take 
2^0 = (2^1 + 2^2) /2. This enables us to do the x-integrals, which yields a couple of (5-functions, 
which we can then integrate over to obtain 



d^k, 



a^ 



((a,a.)B) = {CkXk2)o + ^j -^,mB)a-kB)) ^,[e-'^- {ae-^'k{)ae-^-k2)) ^]\^^^^ 

and using the definition (CfciCfc2) = {27i)^6^{ki + k2)P(;{ki), we obtain the relation 



1 r ^3 



{{CkXk2)B) = (aia.)o + 2(2^)'^'(^i + ^2)/ ^Pc(M(3 + ^j^m(^i; 



d^ks „ ,, X ,„ d 



1 r rl^k 

(a,a.)o + (2 K - 1)' + as){2n)'S%h + ^2)^(^1) J (^n(^^)(B-18) 
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Above, working for simplicity to leading order in slow-roll, we have included the contribution 
from the running of the spectral index Ug = dns/d\og{k) from applying (3 + d/d\og{k)) the 
second time on (n^ — l)P^(/ci). Note that to higher order in slow-roll, other new interesting 
non-scale invariant terms proportional to log{k) begins to appear. 
Now we can use 



(c^(-)>. 

to obtain to leading order in slow-roll 



a.H, (2vr)' 



Mk 



B) 



(B.19) 



(feCg.)^) = {C^,S>^,), + ilin,-ir + ^s)i2nr6'ih + h)Pdh){eix))^ (B.20) 



in agreement with eq (4.12). 



C Loop calculation details 

The steps to calculate the quantities Afe(^), Bkijf), and Ckijf) are outlined in the main text. We 
do not include the intermediate steps as they result in somewhat unwieldy expressions, which 
were calculated via Mathematica. However, a more complete result, expanded in powers of the 
UV and IR cutoffs, is: 



Ak{7]) + Bk{r]) 





H^ H^ J 


"4 2 / 1 \ 101" 
[3log(*)-21og(A„)--log(^^J+g^J 




2A;3 (27r)2 )^ 


+ 


{kvf 


'2 2 2 / 1 \ 241" 
log{k) \og{AiB) ^ log L + ..n 
15 5 3 \Auv J 450 




+ 


{kvr 


\> 


Mi/vA , 4 1 

a- 4- 

"^y k J 75_ 


1 V ^^uv / 





(c.i: 



The quantity Ck can also be computed: 

H^ H^ f dq 



Ckiv) 



3fc4 (27r) 



Q 



.PgY - kY - -kqW 



2k 



3k 



(C.2) 



which must be UV and IR regulated. 

Likewise, we have computed the analogous integrals A'^, 5^, and C[ arising from the corre- 
lator ((Tyfc^ dfcj ) . The sum of the former is 



A'M + B'M 



2fc3(27r)2 
"4 



+ {kr]) 



15 



{kvr 

log 



15 



log(A 



IR) 



^°Ka^ 



7 



4 



k 



A 



UV 



+ 



38 



4k 



225 15A 



UV 



+ A 



450 5 
1 



log(A;) 



>■//?' A 2 



A 



UV 



(C.3) 



28 



and the latter is 



H H j dq / , 5 2„6 I ;^5„4 ^;.3„2 4 



^«"> = sj (2^ y 7 1*'" ^ * " - 2* ' " I ■ <•=■*' 
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